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' ABSTRACT 



o 



We investigate the regularity of cluster pressure profiles with REXCESS, a representative sample of 33 local (z < 0.2) clusters 
drawn from the REFLEX catalogue and observed with XMM-Newton. The sample spans a mass range of 10 I4 M o < M500 < 10 15 M o , 
where M soo is the mass corresponding to a density contrast of 500. We derive an average profile from observations scaled by mass 
and redshift according to the standard self-similar model, and find that the dispersion about the mean is remarkably low, at less 
i-Q ■ than 30 per cent beyond 0.2R 50O , but increases towards the centre. Deviations about the mean are related to both the mass and the 

thermo-dynamical state of the cluster. Morphologically disturbed systems have systematically shallower profiles while cooling core 
systems are more concentrated. The scaled profiles exhibit a residual mass dependence with a slope of ~ 0.12, consistent with that 
expected from the empirically-derived slope of the Msoo-ix relation; however, the departure from standard scaling decreases with 
radius and is consistent with zero at /?5oo- The scatter in the core and departure from self-similar mass scaling is smaller compared 
to that of the entropy profiles, showing that the pressure is the quantity least affected by dynamical history and non-gravitational 
physics. Comparison with scaled data from several state of the art numerical simulations shows good agreement outside the core. 
Combining the observational data in the radial range [0.03-1] ^500 with simulation data in the radial range [l^t]-^5oo, we derive 
a robust measure of the universal pressure profile, that, in an analytical form, defines the physical pressure profile of clusters as a 
' function of mass and redshift up to the cluster 'boundary'. Using this profile and direct spherical integration of the observed pressure 

profiles, we estimate the integrated Compton parameter Y and investigate its scaling with M 500 and L x , the soft band X-ray luminosity. 
£f~) , We consider both the spherically integrated quantity, Y sp h(R), proportional to the gas thermal energy, and the cylindrically integrated 

quantity, Y Ly i(R) = Y SZ D\, which is directly related to the Sunyaev-Zel'dovich (SZ) effect signal. From the low scatter of the observed 
^sph(^50o)-yx relation we show that variations in pressure profile shape do not introduce extra scatter into the K sp h(^5oo)-A^50o relation 
as compared to that from the Y x -M soo relation. The Y sph (R 500 )-M 500 an d ^sph(^50o)-^x relations derived from the data are in excellent 
agreement with those expected from the universal profile. This profile is used to derive the expected Ksz-Msoo and Ysz-Lx relations 
for any aperture. 

Key words. Cosmology: observations, Cosmology: dark matter, Galaxies: cluster: general, (Galaxies) Intergalactic medium, X-rays: 
galaxies: clusters 

x. 

H , 1. Introduction As the gas pressure is directly related to the depth of the 

- - gravitational potential, Ys Z D 2 . is expected to be closely related 

Galaxy clusters provide valuable information on cosmology, t o the mass. Numerical simulations (e.g.. Ida Silva et al. 2 004: 
from the nature of dark energy to the physics driving galaxy and | Mo tl et al. 2005I iNaeai 2006HBonaldi et al. 2007I) and analvti- 
structure formation. Clusters are filled with a hot ionised gas that ca i mode i s dReid & Soereel 2006I) of cluster formation indicate 
can be studied both in X-ray and through the thermal Sunyaev- that the intrinsic scatter of the y sz -M rela tion is low, reg ard- 
Zel'dovich (SZ) effect, a spectral distortion of the cosmic mi- less of the cluster dynamical state (see also IWik et al. 200i or 
crowave background (CMB) generated via inverse Compton the exact details of the gas physics . However, the normalisation 
scattering of CMB photons by the free electrons. Its magnitude f the relation doe s depend on the gas physics dNaeai 2 006: 
is proportional to the Compton parameter y, a measure of the gas | Bo naldi et al. 20071) . as does the exact amount of scatter , the de- 
pressure integrated along the line-of-sight, y = (o- r /m e c 2 ) J Pdl, ta n s f w hi c h are still under debate dShaw et al. 2008I) . Given 
where cr T is the Thomson cross-section, c the speed of light, m e mat this relation, and the underlying pressure profile, are key in- 
the electron rest mass and P = n e T is the product of the electron gre dients for the use of on-going or future SZ cluster surveys for 
number density and temperature. The total SZ signal, integrated cosmology, and provide invaluable information on the physics of 
over the cluster extent, is proportional to the integrated Compton the intra-cluster medium (ICM), it is important to calibrate these 
parameter Y sz , Y SZ D 2 A = (o- T /m e c 2 ) j PdV, where D A is the an- quantities from observations, 
gular distance to the system. 

In recent years, SZ observational capability has 

Send offprint requests to: M. Arnaud, e-mail: made spectacular progress, from the first spatially re- 
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dPointecouteau et al. 199911 200 U iKomatsu et al. 199911 20011 ) 
to the first discovery of new clusters with a blind SZ survey 
dStaniszewski et al. 20091) . Spatially resolved SZE observa- 
tions directly probe the mass weighted temperature along the 
line of sight. By contrast, temperatures derived from X-ray 
spectra, by fitting an isothermal model to a multi-temperature 
plasma emission along the line of sight, are likely to be 
biased (Mat hiesen & Ev rard 2001). Alth ough schemes to 
correct for this effect have been defined ([Mazzotta et al. 2004; 
IVikhlinin 20061) . it remains a potential source of systematics. 

Stacking analysis of WMAP data around known X- 
ray clusters has allowed statistical detection of a scaled 
pressure profile (Af shordi et al. 20071) or a spatially re- 
solved decrem ent (Lieu et al. 2006; At rio-Barandela et al. 2008t 
Diego & Partridge 2009), showing clear discrepancies with the 
prediction of a simple isothermal y6-model. Pressure or temper- 
ature profiles of individual clusters have started to be derived 
from combined analysis of X-ray and SZ E imaging data, us- 
ing non-parametric deprojection methods (Nordet al. 2009) or 
more realistic models th an the j6-model (Kitavama et al. 20041; 
Mro czkowski et al. 20091) . Interestingly, the profiles are found 
to be consiste nt with profiles derived using X-ray sp ectroscopic 
data (see also iJia et al. 2008; lHalverson et al. 20 09). However, 
such studies are still restricted to a few test cases, particularly 
hot clusters. 

The Ysz—M relat ion has been recently derived by 
Bonamen te et al. (20 08). an important step forward as com- 
pared to previous work based on centra l decrement measure- 
ments using heter ogenous data sets dMcCarthv et al. 2 003; 
iMorandi et al. 2 007); however, quantities were estimated within 
^2500 ~ 0.44,^500 El and assuming an isothermal yS-m odel, 
which may provide a biased estimate dHallman et al. 20071) . In 
addition, the first scaling relation using weak lensing masses, 
rather than X-ray hydrostatic masses, has now appeared 
dMarrone et al. 20091) . although constraints from these data are 
currently weak. 

In this context, statistically more precise, albeit indirect, 
information can be obtained from X-ray observations. A key 
physical parameter is Yx, the X-ray analogue of the integ rated 
Compton parameter, introduced by iRravtsov et al. (20061) . Yx 
is defined as the product of M g ,5oo, the gas mass within R500 
and Tx, the spectroscopic temperature outside the core. The lo- 
cal Msqq -Yx relation for r e laxed clusters h a s recently been cal- 
ibrated dNagai et al . 20071 iMaughan 20071; lArnaud et al. 20071; 
Vikhlinin et al. 2009), with excellent agre ement achieved be - 
tween various observations (e.g., see lArnaud et al. 20071) . 
However, the link between Yx and Fsz depends on cluster struc- 
ture through 



<x x 1 



m e c 2 fi e m p (n e ) Rm T x 



(1) 



where the angle brackets denote volume averaged quantities. 
From Eq. Q] it is clear that an understanding of the radial pres- 
sure distribution and its scaling is important not only as a probe 
of the ICM physics, but also for exploitation of these data. High 
resolution measurements of the radial density and temperature 
distribution are now routinely available from X-ray observations 



1 Here and in the following, M e and R s are the total mass and radius 
corresponding to a density contrast, 6, as compared to p c (z), the critical 
density of the universe at the cluster redshift: M s = (4n/3)Sp c (z)Rl- 
M500 corresponds roughly to the virialised portion of clusters, and is 
traditionally used to define the 'total' mass. 



but the pressure profile structure and scaling have been relatively 
little studied. The pressure pro file s of groups have bee n stud- 
ied bv lFinoguenov et aL (2006) an dlJohnson et al. (20091) . In the 
cluster regime, iFinoguenov et al. (2005b analysed the 2D pres- 
sure distribution in a flux-limited sample of 6 hot (kr > 7keV) 
clusters at z ~ 0.3 showing fluctuations at the 30% level around 
the mean profile, scaled by temperature. To our knowledge, 
the only study of pressure profiles scaled by mass is that of 
iNagai et al. (20071) . who used Chandra X-ray observations to 
derive a universal pressure profile, with the external slope de- 
rived from numerical simulations. However, their sample was 
restricted to hot (kT > 5keV) relaxed clusters, which are all cool 
core systems, and contained five objects. For the reasons men- 
tioned above, it is of considerable interest to extend this analysis 
to data from a larger and more representative sample of the clus- 
ter population. 

In this paper we do this by investig ating the regularity o f 
cluster pressure profiles with REXCESS (Bohrin ger et al. 2 007). 
a representative sample of 33 local (z < 0.2) cluste rs drawn 
from the REFLEX catalogue (Bohringer et al. 2004) and ob- 
served with XMM-Newton. We derive an average profile from 
observations scaled by mass and redshift according to the self- 
similar model and relate the deviations about the mean to both 
the mass and the thermo-dynamical state of the cluster (Sec. [3). 
Comparison with data from several state of the art numerical 
simulations (Sec. |4]i shows good agreement outside the central 
regions, which is the most relevant aspect for the Isz estimate. 
Combining the observational data in the radial range [0.3-1] Z?5oo 
with simulation data in the radial range [1-4] 7?5oo allows us to 
derive a robust measure of the universal pressure profile up to the 
cluster 'boundary' (Sec. [5]). Using this profile or direct spheri- 
cal integration of the observed pressure profiles, we estimate the 
spherically and cylindrically integrated Compton parameter and 
investigate its scaling with Yx, M^qq and Lx, the soft band X-ray 
luminosity (SecO. 

We adopt a ACDM cosmology with Hq = 70 km/s/Mpc, 
Q.M - 0.3 and Qa = 0.7. h(z) is the ratio of the Hubble con- 
stant at redshift z to its present value, Hq. Tx is the tempera- 
ture measured in the [0.15-0.75] R500 aperture. All scaling rela- 
tions are derived using the B CES orthogonal regressio n method 
with bootstrap resampling dAkritas & Bershadv 19961) . and un- 
certainties are quoted throughout at the 68 per cent confidence 
level. 

2. The rexcess data set 

A description of the REXCESS sam ple, including XMM-Ne wton 
observation details, can be found in lBohringer et al. (2 007). The 
two clusters RXCJ0956.4-1004 (the Abell 901/902 superclus- 
ter) and J2 157. 4-0747 (a bimodal cluster) are excluded from 
the present analysis. Clus ter subsample cla ssification follows 
the definitions described in lPratt et al. (2009ah : objects with cen- 
tre shift parameter (w) > O.OI.R500 are classified as morpho- 
logically disturbed, and those with central density h(z)~ 2 n e $ > 
4 x 10~ 2 cirT 3 as cool core systems. 

The gas den sity profiles, njj), were derived by 
ICroston et al. (20081) from the surface brightness profiles 
using the non-parametric deprojection and PSF -deconvolution 
technique introduced by ICroston et al. (20061) . The density 
at any radius of interest is estimated by interpolation in the 
log-log plane. The pro cedure to extract the 2D temperature 
profiles is detailed in Prat t et al. (2009al) . The 3D profiles, 
T(r), were derived by fitting convolved parametric models 
d Vikh linin et a l. 20061) to these data, taking into account projec- 
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Fig. 1. The pressure profiles of the REXCESS sample. Pressures are 
estimated at the effective radii of the temperature profile (points with 
errors bars). A line connects the data points for each cluster to guide the 
eye. The data are colour coded (from cold-blue to hot-red) according 
to the spectroscopic temperature, T x . 



tion and PSF effects dPointecouteau et al. 20 04) and weighting 
the contri bution of tempe rature components to each ring as pro- 
posed by IVikhlinin (20061) to correct for the spectroscopic bias 
mentioned above. A Monte Carlo procedure is used to compute 
the errors, which are then corrected to take into account the 
fact that parametric models over-constrain the 3D profile. Full 
details will be given in a forthcoming paper. As the temperature 
profiles are measured on a lower resolution radial grid than 
the density profiles, the pressure profiles, P(r) - n & (r)T(r) , 
are estimated at the weighted effective radii dLewis et al. 2 003) 
of each annular bin of the 2D temperature profiles. They are 
presented in Fig.Q] 

Since the sample contains systems in a variety of dynam- 
ical states, we choose to use Yx as a mass proxy rather than 
the hydrostatic mass. Extensive discussion of how this could 
affect our results is presented in Sec. 13.41 For each cluster, 
M500 is estim ated iteratively from the Mgoo-ix relation, as 
described in iRravtsov et al. (2006I) . We used the updated cal- 
ibration of the Msp p-Fx relation, obtained by combining the 
I Arnaud et al. (2 007) data on nearby relaxed clusters observed 
with XMM-Newton with new REXCESS data (Arnaud et al., 
in prep). The samp le comprises 20 clusters: 8 clusters from 
lArnaud et al. (20071) . excluding the two lowest mass clusters 
whose M500 estimate requires extrapolation, and the 12 relaxed 
REXCESS clusters with mass profiles measured at least down to 
6 = 550. The derived Msoo-Fx relation 



h(z) 2 ' 5 M : 



10 1 ' 



2x 10 14 h;„ 5/2 M G keV 



h™ 1 M Q (2) 



is consistent with the relation derived bv l Arnaud et al. (2007) but 
with improved accuracy on slope and normalization. 

The slope differs from that expected in the standard self- 
similar model (a = 3/5) by only ~ 2<x. We will thus also con- 
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Fig. 2. The scaled pressure profiles of the REXCESS sample, colour 
coded according to the (thermo)dynamical state (see labels and Sec. IS. 
The radii are scaled to R soo and the pressure to P 50O as defined in Eq.[5] 
with M500 estimated from the Msoo-Yx relation (Eq.O. Full lines: pres- 
sure profile as in Fig. Q] with data points omitted for clarity. Dotted 
lines: extrapolated pressure (see text). The thick black line is the aver- 
age scaled profile and the grey area corresponds to the ±l<r dispersion 
around it. Middle panel: ratio of the average profile of cool core (blue) 
and disturbed (red) systems to the overall average profile. Bottom panel: 
The solid line is the statistical dispersion as a function of scaled radius. 
Dotted line: additional dispersion expected from the intrinsic dispersion 
in the A/500-Ix relation. Dash-dotted line: quadratic sum of the two dis- 
persions. Dashed line: dispersion obtained for M500 estimated from the 
standard slope Msoo-^x relation (Eq.[3]l. 



sider the A/500-Fx relation obtained by fixing the slope to its 
standard value: 



h(zf 15 M 500 = 10 1 



2 x 10» 4 h TO 5/2 M G keV 



3/5 



3. Scaled pressure profiles 

3.1. Scaled profiles 
The scaled pressure profiles 
P(r) 



p(x) 



where x = 

P500 R500 



(4) 



are presented in Fig. [2] The pressure is normalised to the char- 
acteristic pressure P500, reflecting the mass variation expected 




Fig. 3. The scaled density (left panel) and temperature (middle panel) profiles of the REXCESS sample. Each profile is colour coded according 
to the cluster (thermo)dynamical state (see labels and SecO. The radii are scaled to Rsqo, estimated from the M 500 -Y x relation (Eq.[2jl. The density 
is scaled to the mean density within R$oo and the temperature to Tx, the spectroscopic temperature measured in the [0.15-0.75] R500 aperture. In 
each panel, the thick black line is the average scaled profile. Right panel: logarithmic deviation of the scaled temperatures from the average scaled 
profile versus the corresponding deviation for the density, at each effective radius of the temperature profile annular bins. Data corresponding to 
r/Rsoo < 0.2 and r/R^oo > 0.2 are marked with filled and open circles, respectively. The deviations are anti-correlated in the core. 



in the standard se lf-similar model, purely based on gravitation 
dNagai et al. 20071 and Appendix©. 



500 



1.65xl0~ 3 /i(z) 



8/3 



M< 



500 
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14U-1 



h^o M 



2/3 



h^keV cirT 3 (5) 



For comparison we also plot in Fig. [3] the scaled temperature 
profiles, t(x) = T(r)/Tx as well as the scaled density profiles, 
n e (x). Note that the density profiles have been normalised to the 
mean density within /?5oo, so that the dispersion is only due to 
variations in shape Q. 

The resolution in the centre and radial extent of the pressure 
profiles are determined by that of the temperature profiles, in 
practice the effective radius of the inner and outer annular tem- 
perature profile bins, which varies from cluster to cluster (see 
Fig. |2|. In particular, the peaked emission of cool core clusters 
allows us to measure the profiles deeper into the core than for 
disturbed clusters, which have more diffuse emission (see also 
Sec.|33]). 

3.2. Average scaled pressure profile 

We computed an average scaled pressure profile, p(x), from the 
median value of the scaled pressure in the radial range where 
data are available for at least 15 clusters without extrapolation 
(about [0.03-1] /?5oo)- However, to avoid a biased estimate of 
the average profile in the core, where the dispersion is large and 
more peaked clusters are measured to lower radii (Fig. [2), it is 
important to include all clusters in the computation. For this pur- 
pose, we extrapolated the pressure profiles in the core using the 
best fitting temperature model used in the deprojection of the 
temperature profile. This extrapolation is only weakly model de- 
pendent since it essentially concerns disturbed clusters (Fig. [2]i 
which are observed to have rather flat central temperature pro- 
files (Fig. [3). The average profile is plotted as a thick line in 
Fig. |2] The dispersion around it is defined as the plus or minus 
standard deviation from the average profile, computed in the log- 
log plane. 



2 The normalisation of the density profiles, scaled according to 
the standard self-s imilar model, varies with mass as shown by 



ICroston et arT2 008). 



3.3. Dispersion, radial structure and dynamical state 

For a perfectly self-similar cluster population, the scaled pro- 
files should coincide. The dispersion around the average scaled 
profile is less than 25% beyond the core (r > 0.27?soo) and in- 
creases towards the centre (bottom panel of Fig. [2j. This disper- 
sion reflects a variation of shape with cluster (thermo)dynamical 
state, as clearly seen in Fig. [2] shallower profiles, at all radii, 
are observed for morphologically disturbed clusters while the 
cooling core clusters have the most concentrated profiles. The 
typical difference between the average profiles of these two pop- 
ulations is ~ 20% in the outskirts and as high as a factor of four 
at 0.03 R 500 (Fig.|2]middle panel). 

When compared to the density profiles (Fig. [3] left panel) 
the pressure profiles are distinctly more regular and present less 
dispersion in the core. The reason lies in the anti-correlation be- 
tween the deviation of scaled temperatures and densities from 
their respective average scaled profiles, m(x) and i(x), as shown 
Fig. [3](right panel). For data interior to r < 0.2/?5oo, a Spearman 
rank test finds a probability of 10~ 7 that the anti-correlation be- 
tween log(n^(jc) - \og{x\(x)) and log(f(x) - log(lt(x)) occurs by 
chance. The correlation disappears at large radii (probability of 
0.6 for r > 0.2/?5oo). Qualitatively, this is the result of the well- 
known fact that cool core clusters have peaked density profiles, 
with a temperature drop in the centre, while unrelaxed objects 
have flatter density cores (Fig. [3] left panel) and constant or in- 
creasing temperature toward the centre (Fig. [3] middle panel). 



3.4. Dependence on mass and mass-proxy relation 

Since we derived M500 from the Msoo-ix relation, the scaling 
quantities R^oq and P500 and the pressure profiles are not in- 
dependent, as they are both related to the product of the gas 
density and temperature. We first examine how this may affect 

our results. From the definition of the pressure P(r)=n e (r)T (r), 

2/3 

and noting that P500 °c MJ m and that Yx — Mg ^ooTx °c 
(rieir^R^R^Tx K {n e (r)) R5m T x M 5 o( h where the angle brack- 
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Fig. 4. The scaled pressure profiles of the REXCESS sample, colour coded according to spectroscopic temperature measured in the [0.15-0.75] 
Rioo aperture (left panels). Right panels: corresponding scaled pressure estimated at different values of scaled radii plotted as a function of cluster 
mass. Full lines: power law fit at each scaled radius, with the best fitting slope given in the labels. Dotted line: value for the average scaled profile 
at that radius. The cluster masses, M500. are estimated from the Msoo-Yx relation, either the best fitting empirical relation (Eq.(2] top panels) or the 
relation obtained from fixing the slope to its standard value (Eq.[3] bottom panels). 



ets denote a volume average within Rsoo, the scaled pressure 
p(x) = P(xRsoo)/Psoo is proportional to 

, . P(xR 500 ) (ne(r)T(r)) Rm Y x 

This equation makes explicit the link between the scaled pres- 
sure profiles and the Msoo-Jx relation. The first two dimension- 
less terms in the right hand part of the equation purely depend 
on the internal gas structure within R500. They determine the av- 
erage shape of the scaled profile. The third term depends on the 
global cluster scaling properties between Yx and M500 and de- 
termine both the normalisation of the average scaled profile and 
the 'typical' mass dependence of the profiles (discussed at the 
end of the section). 



Using M500 values derived from the M^qq-Yx relation, rather 
than the 'true' M500 value, is equivalent to assuming a perfect 
correlation between M500 and Fx, i-e with no scatter. Provided 
that the correct A/500-lx relation is used and that o~i 0S) my does 
not depend on mass or dynamical state, use of the M 5 qq-Y x re- 
lation will not introduce a systematic bias into the scaled pro- 
files, but their dispersion will be underestimated. Let us de- 
fine the intrinsic scatter of the A/500-lx relation, o~i 0S) my> as 
the standard deviation of log(M5oo) from the value from the 
best fitting relation at a given Yx- We can estimate the addi- 
tional dispersion due to o~ i g,MF from the effect on the average 
scaled profile of a variation of log(M5oo) by ±<t\ ^my- Since 

R$oo oc Mjqq and P500 K ^500 > me profile is translated in 
the log-log plane by ±l/3cH og jW y and +2/3cr log M y along the x 
and y axis, respectively. Assuming cr\ os jay = 0.04 (about 10%, 
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iKravtsov et al. 2006t lArnaud et al. 20071) . the additional disper- 
sion (in dex units), computed from the difference between the 
translated profiles at a given scaled radius, is plotted in the bot- 
tom panel of Fig. [2] It is non-negligible beyond the core, but the 
total dispersion, estimated by summing quadratically this addi- 
tional contribution, is expected to remain below 30%. It is negli- 
gible in the core, where the dispersion is dominated by structural 
variations. 

Finally, the Msqq-Yx relation being derived from mass es- 
timated using the hydrostatic equilibrium, we expect an offset 
between that relation and the 'true' Msm-Yx relation. The M500 
used in this study are thus likely to be underestimated. The effect 
of such a bias is to translate all the scaled profiles together (pro- 
vided that it is a simple factor independent of mass). This will 
not affect any shape or dispersion analysis but the normalisation 
of the mean scaled profile will be biased high. This is further 
discussed in Sec. l4.3l and Sec. 16.21 

We now turn to the question of the variation of the pres- 
sure profile normalisation with mass. From the definition of 
P500, any deviation from the standard self-similar scaling will 
appear as a variation of the scaled profiles p(x) with mass, 
p(x) = p(x, M500). It will also translate into a non-standard slope 
Q'MYx f° r the M500— Jx relation. From Eq.[6]we expect that p(x) 
increases slightly with mass as Yx/M^, i.e as M^ jQ with: 



a P = 



ffMYx 



0.12 



(7) 



for the best fitting slope ayiY x - 0.561 (Eq. [2J. We show in 
the left- top panel of Fig. H] the scaled profiles colour coded as a 
function of Tx- There is some indication that hotter (thus more 
massive) clusters lie above cooler systems. To better quantify 
the variation with mass, the right-top panel of the figure shows 
the variation with M500 of the scaled pressure, p(x), for different 
scaled radii, x - r/R^oo. At each radius, we fitted the data with a 
power law p(x) oc . At all radii, the slope a(x) is consistent 
with the expected 0.12 value (Fig. HI right top panel), and the 
mean slope is 0.10 ± 0.02. The pivot of the power law, where 
the pressure equals the average scaled value, p(x) = p(x), is 
about M500 ~ 3 x 1O 14 M . In a first approximation, the mass 
dependence of the scaled profiles can then be modelled by: 



p(x, M 500 ) = p(x) 



500 



3 x 10 



141,-1 



h- ( ; m 



«p=0.12 



(8) 



where p(x) is the average scaled profile derived in Sec. 13.21 

However, there is some indication that the mass dependence 
of the profiles is actually more subtle than a global normalisa- 
tion variation. The variation in slope of the power law fits shows 
that the mass dependence decreases with radius, with a(x) = 
0.22 ± 0.16 at r = 0.1R 500 and a(x) = -0.01 + 0.16, consistent 
with zero at /?5oo- In other words, the departure from standard 
scaling, likely to be due to the effects of non-gravitational pro- 
cesses, becomes less pronounced as we move towards the clus- 
ter ou tskirts, behaviour that was also notic ed in the entropy pro- 
files dNagai et al. 20071: iPratt et al. 2009bl) . Note, however, that 
the mass dependence is weaker for the pressure than for the en- 
tropy: the pressure slopes ar e about two time s smaller than those 
of the entropy (Fig. g] and IPratt et al. 2009U their Fig 3). The 
variation of a(x) with x can be adequately represented by the 
analytical expression, a(x) = ap + a' p (x), with: 



(x/0.5) 3 

< W = o.io- ( . P + o.io )T -i^ 



(9) 




0.01 



0.10 
Radius (R 500 ) 



1.00 



Fig. 5. The scaled pressure profiles (green lines) derived from 
IVikhlinin et al. (2 006) Chandra data on relaxed clusters compared to 
the scaled profiles of the R EXCESS sample excluding morphologi- 
cally disturbed clusters (same colour code as in Fig.|2]l. The thick green 
dotted line is the average Chandra profile. Bottom panel: ratio of that 
average Chandra profile to that of REXCESS for all morphologically 
undisturbed objects (dotted line) or only cool core clusters (full line). 



yielding to the more accurate, although more complex, model 
for the scaled profiles: 



p(x, M500) = p(x) 



M 50 



3 x 10 



h^„ M 



(10) 



We then compared to the results obtained using M500 derived 
from the self-similar Msqq-Yx relation with slope 3/5 (Eqj3). The 
scaled profiles are plotted in the bottom panel of Fig. g] In this 
case, we do not expect any dependence of p(x) with M500, and 
this is indeed the case: the slopes ap(x) are consistent with zero 
at all radii (right bottom panel). The dispersion in scaled pro- 
files is also smaller (see Fig. |2]bottom panel). In that case, the 
dispersion is only due to structural variations, while in the non- 
standard case, the mass dependence of p(x) also contributes to 
the dispersion. 

3.5. Comparison to Chandra results for relaxed clusters 

In Figure H] we plot the pressure profiles presented in 
Nagai et al. (2007), derived from Chandra data analyzed by 
Vikhlinin et al. (20061) . We only consider clusters with measured 
M 50 o values, excluding MKW4 (T = 1.4keV) and A2390 (z = 
0.23) which fall outside the Tx and z range of REXCESS, re- 
spectively. We used the published M500 values, derived from the 
hydrostatic equilibrium (HSE) equation, and computed the pres- 
sure from the best fitting par ametric models of th e density and 
temperature profile given in Vikh linin et al. (20"06l) . in the radial 
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range of the observed temperature profile. Since the Chandra 
data set only contains relaxed clusters, they are compared to the 
REXCESS profiles excluding morphologically disturbed objects. 

All Chandra profiles, except one, lie within the range of the 
REXCESS profiles. The bottom panel of Fig. |5]shows the ratio 
of the average Chandra profile to the average REXCESS pro- 
file. The agreement between the average profiles, both in shape 
and normalisation, is nearly perfect beyond the core, where the 
dispersion of the scaled profiles is lower. However, on average, 
the Chandra profiles are slightly more peaked towards the cen- 
tre (dotted line in bottom panel of Fig. (5j and have a smaller 
dispersion than the 'relaxed' REXCESS clusters. Better agree- 
ment is found with the average REXCESS profile for cool core 
clusters (full line in bottom panel of Fig. |5j. This is not surpris- 
ing, since all clusters in the Chandra data set present the central 
temperature drop characteristic of cool core clusters. 

This good agreement is an indication of the robustness of 
scaled pressure profile measurements with current X-ray satel- 
lites. The comparison also illustrates the importance of consider- 
ing a representative cluster sample to measure the average profile 
and dispersion in the core. 



4. Comparison with numerical simulations 

4.1. The data set 

We consider three large samples of simulated clusters 
at redshift zero extracted from ACDM cosmological N- 
body/hydrodynamical simulations (Qm = 0.3, Qa - 0.7). The 
data set includes the samples from Borgani et al. ( 2004. 
hereafter BO), iPiffaretti & Valdarnini (2008[ PV) and 
Naga i et al. (2007I NA). All simulations include treatment 
of radiative cooling, star formation, and energy feedback from 
supernova explosions. The three simulated data sets are fully 
independent and derived using different numerical schemes and 
implementations of the gas physics (see references above for 
full description). This allows us to check the robustness of the 
theoretical predictions of the pressure profiles by comparing the 
three simulated data sets. The fact that the NA simulation was 
undertaken on a mesh-based Eulerian code, while the PV and 
BO simulations were derived from particle-based Lagrangian 
codes is particularly relevant, considering some well known 
cluster-scale discrepancies between the numeri cal approaches , 
such as is seen in the entropy profiles (see, e.g., IVoit et al. 2005t 
Mitc hell et al. 2 009, and references therein). The star formation 
algorithm and the SN feedback model are also quite different 
both in implementation and in feedback efficiency. 

In order to avoid comparison with inappropriately low mass 
objects we impose the REXCESS lower mass limit M500 > 
1O 14 M , leading to a final number of simulated clusters of 93, 
88, and 14 for the BO, PV, and NA samples, respectively. We 
computed the pressure profile for each cluster using the mass- 
weighted gas temperature, since the deprojection of the observed 
profile takes into account the spectroscopic bias (Sec. The 
assumed baryon densities are Qb = 0.039, 0.049, 0.043 for the 
BO, PV, and NA samples, respectively. The assumed baryon 
fraction, fb = Q^/Q,,, has a direct impact on the gas density 
and thus pressure profile at a given total mass. We thus cor- 
rected the gas profiles by the ratio between th e assumed fb value 
and the WMAP5 value dDunklev et al. 2009h for each sample. 
To scale each individual pressure profile we consider both the 
'true' 



DSim 

MHSE _ il^HSE/ - pHSE\ 
500 - M (< ^500 > 



and M^™ values and the hydrostatic values and 
The former are derived from the total 
mass distribution in the simulation. The latter was derived from 
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Fig. 6. The scaled pressure profiles derived from num erical simu- 
lations of iBorgani et al. (20 04) (pink), iNagai et al. (2007f) (blue) and 
PifFaretti & Valdarnini (2008) (green). Black line: overall average pro- 
file (see text), coloured lines: average profile for each simulation with 
the coloured area corresponding to the dispersion around it. Bottom 
panel: ratio of each simulation average profile to the overall average 
profile. 



the gas density and temperature profiles and the hydrostatic equi- 
librium equation, using the same procedure for all c lusters. As 
in previous work (e.g., Piffaretti & Valdarnini 2008, and refer- 
ences therein), we find that underestimates the true mass. 
We find a mean bias for the whole sample of - 1 3 per cent with a 
dispersion of + 1 6 per cent; the average bias estimated for the dif- 
ferent simulations agrees within a few percent at all radii larger 
than 0.1/? 5 no- 



4.2. Comparison of numerical simulations 

We derive the average scaled profile for each simulation, and 
the dispersion around it, from the median value and 16 and 84 
percentiles of the scaled pressure distribution at a given scaled 
radius. We also compute an average simulation profile. Since the 
average profile computed from the total sample would be biased 
by the number of objects in the largest data set, we average the 
three mean profiles from each simulation data set, and calculate 
the dispersion from all available profiles. The results derived us- 
ing the true mass are shown in Fig. [6] 

Taking into account that the profiles vary by more than 5 or- 
ders of magnitude from the cluster centre to the outskirts, the 
agreement between the three simulations is exceptionally good. 
The profiles agree within 20% between ~ 0.1 and ~ 3 /?™ (Fig. 
|6]lower panel). As expected, larger differences are found in the 
core, where non-gravitational processes are more important and 
where the differences in their implementation in the codes will 
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Fig. 7. Comparison of the R EXCESS scaled profiles with the pre- 
diction of numerical simulations. Black lines: R EXCESS data (as in 
Fig.[2]l. Thick black line: average REXCESS scaled profile. Red line: 
average simulation profile and dispersion around it (orange area) using 
the hydrostatic mass. Dotted red line: same using the true mass. Bottom 
panel: ratio of these average simulation profiles to the REXCESS av- 
erage profile. 



only relaxed clusters, the median bias on M500 changes by only 
2%, the main effect being a factor of 2 decrease in its dispersion. 

The simulation prediction and the REXCESS data agree well 
in the external part (r > Q.2R 5 qo), with the observed profiles ly- 
ing within the dispersion around the average simulation profile 
(Fig. |7]). Remarkably, the observed and simulated average pro- 
files are parallel above O.4R500 (i.e they have the same shape), 
with a normalisation offset of only ~ 10% (Fig. [7] bottom panel). 
The slight underestimate of the pressure in the simulations is 
similar to the offset observed for the M^q-Yx relation and may 
be due, at least in part, to over-co ndensation of hot gas in the 
cold dense phase (see discussion in lArnaud et al. 20O7I) . As we 
move towards the centre, the agreement progressively degrades, 
the simulations predicting more peaked profiles than those ob- 
ser ved (Fig. [7] bottom panel). This behaviour was also noticed 
by iNagai et al. (20071) when comparing their simulations with 
Chandra relaxed clu sters, and it is a lso observed for the tem- 
perature profiles (see lPratt et al. 20071) . As mentioned above, the 
core properties are most sensitive to non-gravitational processes 
and these discrepancies are again likely to reflect the fact that 
modelling of the processes is still inadequate. 

The average simulation profile derived using the true mass 
for each simulated cluster is also shown in the figure (dotted 
lines). As compared to the scaling based on R^^ and M^ E , 
the scaled profile of each cluster is translated to the left and to 
the bottom in the log-log plane. The average profile lies below 
the profile based on the hydrostatic values, as expected from the 
mean bias between and M^^. The offset with the observed 
profile in the outer region becomes more significant, about 30%. 

In conclusion, there is an excellent agreement in shape be- 
tween the simulated and observed profiles for the cluster outer 
regions, which is the most relevant aspect for the Fsz estimate. 
The better agreement in normalisation with the simulations when 
using the hydrostatic mass suggests that the hydrostatic X-ray 
masses used to scale the observed profiles are indeed underesti- 
mated. 



become more evident. The BO profiles are available only up to 
the 'virial' radius, 2.03 Rf™ but the PV and NA profiles are 
traced up to 10/?^™, where they deviate significantly, but still 
agree within the dispersion. However, the differences are sytem- 
atic with the PV profiles lying below the NA profiles. This may 
hint at a difference in the way in which Lagrangian and Eulerian 
codes behave in the IGM-WHIM regime. Note also the flatten- 
ing of the pressure profile in the outskirts, around 5 Rf^, which 
is likely to define the actual boundary of the cluster, where it 
meets the intergalactic medium. In the following we will use this 
boundary to compute the total integrated SZ signal, Fsz- In spite 
of the difference in the pressure in the outskirt, there is good 
agreement on Fsz between the simulations: the SZ signal within 
5 /?5oo computed from the average PV and NA profiles differ by 
-15%, and +9%, respectively, from the value computed using 
the average simulation profile. 

4.3. Comparison of rexcess profiles with simulations 

Figure [7] compares the observed scaled profiles with the predic- 
tion of the simulations. We first consider the simulated profiles 
scaled using the hydrostatic quantities and M^ } y E , since the 
observations rely on hydrostatic mass estimates. Note that we 
used the M500— ix relation calibrated from a sample of relaxed 
clusters, while for the simulations we used ^ S E and M^ E for 
the whole sample. However we checked that, when considering 



5. The universal pressure profile 

As pointed out by Nag ai et al. (20071) . an analytic cluster pres- 
sure profile model is useful both for analysis of SZ observations 
and for theoretical studies. Of prime interest is a model for the 
average scaled profile of the entire cluster population. For nearby 
clusters it can be derived from the present data, the REXCESS 
sample being a representative sample. 

We c onsidered the gene ralized NFW (GNFW) model pro- 
posed bv lNagai et al. (2007b : 

pW ~ 7Z — rr (11) 

(c 5 00^[l+(C500^]^ r)/ff 

The parameters (y,a,/3) are respectively the central slope (r <K 
r s ), intermediate slope (r ~ r s ) and outer slope (r » r s ), where 
r s = ^50o/c5oo, and they are highly correlated with r s . In order 
to constrain the parameters, it is essential to consider a wide ra- 
dial range, including both the core (r < Q.IR^qq) and the cluster 
periphery (r > /?5oo). In particular, /3 remains essentially uncon- 
strained when considering only data within r < R500, resulting 
in large uncertainties in the profile model beyond R 5 qq and thus 
on the corresponding integrated SZ signal. 

Taking advantage of the good agreement between observa- 
tions and simulations in the outer cluster regions, we thus de- 
fined an hybrid average profile, combining the profiles from ob- 
servations and simulations. It is defined by the observed average 
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Fig. 8. GNFW model of the universal pressure profile (green line). 
It is derived by fitting the observed average scaled profile in the ra- 
dial range [0.03-l]/?5oo> combined with the average simulation profile 
beyond R500 (red line). Black lines: REXCESS profiles. Orange area: 
dispersion around the average simulation profile. 



scaled profile in the radial range [0.03-1 ]/?5oo derived in Sec. 13. 21 
and the average simulation profile in the [1— 4]i?5oo region. For 
the simulations, we used the profile based on the hydrostatic 
quantities and renormalised it by +10% to correct for the ob- 
served offset with the observations at r > OAR 500- We fitted 
this hybrid profile with the GNFW model in the log-log plane, 
weighting the 'data' points according to the dispersion. The best 
fitting model is plotted in Fig. [8] with parameters: 

[Po,C50o,y,a,P\ = [8.403 h TO 3/2 , 1.177,0.3081, 1.0510,5.4905] (12) 

Using the dimensionless 'universal' profile, p(x) (Eq. QT| and 
Eq. [T2l . and taking into account the mass dependence estab- 
lished in Sec. 13.41 we can describe the physical pressure profile 
of clusters as a function of mass and redshift (assuming standard 
evolution): 



P(r) = P500 



M 500 



3xlO 14 h^M 



1.65 x 10 _3 ft(z) 



8/3 



M* 



PW (13) 

■i2/3+ap+<4(jt) 



3 x lO^h^Mo 



X |p (x) h 2 ke V cm 3 

with x = r/Rsoo, ap and a' p (x) from Eq. [7] and Eq. [9] and p(x) 
from Eq.Q~T]with parameters from Eq.Q~2] The second term in the 
mass exponent, a P , corresponds to a modification of the standard 
self-similarity (i.e., the steeper mass dependence of the profile), 
while the third term, a' p (x) (Eq. |9), introduces a break in self- 
similarity (i.e., a mass dependence of the shape). The latter is a 
second order effect, which can be neglected in first aproximation. 

We also fitted each individual observed cluster profile with 
the GNFW model, fixing the /? value to that derived above 
(Eq. 112). The best fitting parameters are listed in Appendix ICl 
where we also provide plots of each individual cluster profile 
with its best fitting model. 



6. Integrated Compton parameter scaling relations 

6.1. Definitions and method 

In this section we discuss scaling relations directly relevant for 
SZE studies. We will consider the volume integrated Compton 
parameter Y, for both cylindrical and spherical volumes of in- 
tegration. The spherically integrated quantity, Y sp b_(R), propor- 
tional to the gas thermal energy, is defined as: 



Y sp h(R) 



CTj 



f 

Jo 



AnP(r)r 2 dr 



(14) 



and the cylindrically integrated quantity, Y cy \(R) = Y$zD\, di- 
rectly related to the SZ signal within an aperture 9 = R/Da, is: 

~* b 2 P(r')r'dr' 



(Tj 
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Jo 



2nrdr 



- T sp h(^b) - 



0"t 



r 



(15) 



4nP(r)y/r 2 - R 2 rdr 



where R\, is the cluster radial extent. In the following, we adopt 
Rb = 5/?5oo, as suggested by numerical simulations (Sec. 14.2b . 
Note that the total SZ signal is then equivalently T sp h(5/?5oo) or 

T cy l(57?500)- 

For each cluster, the spherically integrated Compton param- 
eter can be readily computed from the observed pressure profile. 
The F sp h scaling relations can then be directly derived from the 
data for integration radii up to R500, the observed radial range. 
They are presented below in Sec. 16.21 Such a derivation is not 
possible for Y cy \ (or the total Tsz signal): it involves integration 
along the line of sight up to R\, = 5Rsqq, i.e., beyond the observed 
radial range. However, using the universal pressure profile, we 
can compute the volume integrated Compton parameter, Y, for 
any region of interest, and derive the corresponding scaling re- 
lations (presented below in Sec. 16.31 . The two approaches give 
fully consistent results, as shown below. 

Finally, for convenience, we also define a characteristic 
Compton parameter, F500, corresponding to the characteristic 
pressure P500 (see AppendixlAt: 



<X T An 



— R< m P500 



(16) 



2.925 x IQ- 5 'h(z) 
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M 50 o 



3 x 10 



141,-1 
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n 70 Mpc 2 



6.2. Observed Y sph -Y x and y sp h-A^50o relations 

The values for T S ph(^250o) and F sp h(/?5oo), derived from the ob- 
served pressure profiles, are given in Table IC.ll /?2500 is de- 
fi ned as j?2500 = 0-4 4^5qq from the scaling relations presented 
in lArnaud et al. (20051) . The integration was performed using the 
MC deconvolved density and model temperature profiles, allow- 
ing us to propagate the statistical errors, including that on Rsoo- 
We checked that using instead the best fitting GNFW model for 
each profile gives consistent results within the statistical errors. 
Note that the errors on M 5 oo take into account the statistical er- 
rors on the relevant X-ray data, but not the uncertainties on the 
relation itself. The latter are therefore not included in 
the statistical errors on the slope and normalisation of the rela- 
tions. 

Figure [9] shows the T sp h-Tx relations with Fx = Mg,50o7x> 
together with the best fitting power law. We normalised Yx by : 

,-19 M PC 2 
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Fig. 9. The K S p h -Kx relations from REXCESS data. Y sph (R) is the 
spherically integrated Compton parameter, within R2500 (squares) and 
/?50o (circles). Yx = M gi 5oo7x is the product of the gas mass within ^500 
and the spectroscopic temperature Tx- Data points are colour-coded ac- 
cording to cluster dynamical state. Lines: best fitting power law. 



for/i e = 1.148, the mean molecular weight of electrons for a 0.3 
solar abundance plasma. Note that the F sp h-lx relation depends 
only weakly on the assumed Msoq-Yx relation, via the estimate 
of /?5oo only. For some clusters, the computation of F S ph(^50o) re- 
quires extrapolation: by more than 20% for 8 clusters and, in the 
worst case, RXC J2157. 4-0747, the profile of which is measured 
only up to R^t ~ 0.6i?soo, I'sphC^soo) is larger by a factor 1.8 than 
the value within /?det- However, the best fitting 5 / S ph(^50o)-ix re- 
lation is stable to the inclusion or exclusion of clusters requiring 
extrapolation, the best fitting parameters being consistent within 
the errors. 

As mentioned in the introduction, the F sp h-lx relation de- 
pends on the internal cluster structure (Eq. [T). For F sp h(^250oX 
we obtained: 



Isph(^2500) - 10 



-0.272±0.097 



Cxsz Yx 



h -5/2 
"70 



Mpc 2 



1.036±0.020 



h^ 2 Mpc 2 (18) 



The best fitting slope is slightly greater than one (a 2cr effect), 
reflecting the stronger mass dependence of the pressure profile 
in the centre (r < /?250o) as compared the expectation from the 
Msoo-Yx relation (Fig. |4] and Sec. I3.4l i. The intrinsic dispersion 
is Cloggy = 0.054 ± 0.006, with the morphologically disturbed 
clusters lying below the mean relation and the relaxed clusters 
lying above it, a consequence of the shallower profile of the 
former as compared to the latter (Fig. |2). When we move to 
F sp h(^50oX the best fitting slope (1 .003 ± 0.008) becomes consis- 
tent with unity, i.e the shape variation with mass, when averaged 
within 7?5oo, has essentially no effect (see also below). The intrin- 
sic dispersion is no longer measurable, the dispersion is consis- 
tent with that expected from the statistical errors. This is a direct 
consequence of the high similarity of the pressure profiles be- 
yond the core (r > 0.2/?5oo), while the core typically contributes 



o 

CL 



1 



>- 



1 1 1 — 1 — 



- 1 1 1 1 1 — 1 — 



Ycylfeoo) 
'sphA^500/ 

Totol 




Mroczkowski et al., 09 □ 
1 1 1 ■ ■ ■ 1 



10 



15 



M 500 ( h 70 M ©) 



Fig. 10. The Ysz-Msoq relations. Thin black line: power law y sp h(/?5oo)- 
M500 relation best fitting the REXCESS data (black points). Lines: 
y sph (K 5 oo)-Af5oo relation (dashed line), Fcyi(i?5oo) _ -^500 relation (rela- 
tion between Ksz^i within a R500 aperture and M500; full line) and 
Y cyl (5^500 )— M500 relation (relation between the total Y SZ D 2 A signal and 
M500; dash-dotted line) derived from the universal GNFW scaled pres- 
sure profile and for different M500 scaling: standard self-similar scal- 
ing (a = 5/3; green), modified scaling taking into account the non- 
standard slope of the Msoo-lx relation (a = 1.78; blue), and fur- 
ther taking into account the break of self-similarity of the pressure 
profile shape (see text, red). Blue squares: (Foyi(i?5oo). M500) measure- 
ments for 3 clusters: from top to bottom, A1835 (z = 0.25), A1914 
(z = 0.17), and CL J1226 .9+3332 (z = 0.89). They were derived by 
I Mroczkowski et al. (20 09) from a joint analysis of SZA and X-ray ob- 
servations using a GNFW model. 



by less than 10% to F sph (500) (see below and Fig. [TTJ. Fixing 
the slope to one, the best fitting normalisation gives: 



F sp h(^500) 

Cxsz Yx 



0.924 ± 0.004 



(19) 



Note that this ratio is nothing more than the ratio, T mg /Tx, of 
the gas mass weighted t emperature to Tx- It i s less than unity, 
as found in other studies ( Vikhlini n et al. 20"06T) . and as expected 
for decreasing temperature profiles. 

Figure [10] shows the F S ph(^50o)-^500 data together with the 
best fitting relation: 



fc(zr 2/ %h(tf5oo) = 10- 



M,, 



3 x 10 14 h- ' M G 



' Mpc (20) 



Since M500 is derived from the M500— Yx relation, this expression 
does not contain more information than the F sp h(^50o)-F x rela- 
tion, combined with the calibration of the Msoo-Yx relation. As 
expected, the normalisation and slope are consistent with that 
obtained by combining Eq. [2] and Eq [19] and, similar to the 
F sp h(^50o)-F x relation, the scatter is consistent with the statis- 
tical scatter. 

With such a study, based on a mass proxy, we cannot as- 
sess the intrinsic scatter of the 'true' F S ph(^50o)-^500 relation. 
However, we emphasize that our study does show, from the low 
scatter of the F sp h(/?5oo)-ix, that variations in pressure profile 
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shapes do not introduce an extra scatter into the 5 / S ph(^50o)-^50o 
relation as compared to that of the Fx-^500 relation. Since the 
latter was derived from hydrostatic mass estimates using relaxed 
objects, the above 5 / S ph(^50o)-^500 is expected to differ from the 
'true' 5 / S ph(^50o)-A^500 by the offset between the 'true' mass and 
the hydrostatic mass for relaxed objects. 

6.3. Scaling relations from the universal pressure profile 

6.3.1. F sp h-M 5 oo and Y SZ D 2 A -M 500 relations 

Let us first consider F sp h derived from the universal pressure pro- 
file. Combining Eq. El 03 and ED 



Y sp h(R) - *500 



M 500 



3 x lO^h^Mo 



(21) 



f 3/(h, 
Jo 



M500) p(«) u du 



with f(u,M 500 ) = (M500/3 x lO'^Mg)^"'. This term in 
the integral reflects the break of self-similarity in the pressure 
profile (Sec. 0. Neglecting this effect, the corresponding F sp h- 
M500 relation, for any integration radius, is a power law of slope 
a = 5/3 + ffp = 1/ffMYx (Eq.rTlandfTSI). Taking into account this 
effect, the relation is no longer a simple power law. Following 
the behavior of the pressure profiles - a' p (u) decreases with 
radius or equivalently the departure from standard mass scal- 
ing becomes less and less pronounced as we move towards the 
cluster outskirts - the relation is expected to become shallower 
with increasing integration radius, closer to the standard self- 
similar relation (a = 5/3). The relations for various mass scal- 
ings can be compared in Fig. [TO] for F sp h(^50o) and F sp h(57?5oo) 
(i.e the total Fsz signal). The effect of the self-similarity break 
is small. In the mass range M500 = [1O 14 M , 1O 15 M ], F sph 
varies, as compared to the value computed neglecting this effect, 
by [-7%, +8%],[-l%, +0.5%] and [+6%, ~6%], for an integra- 
tion radius of /?2500, ^500 and 5/?soo, respectively. When taking 
into account the self-similarity break, the corresponding effec- 
tive slopes of the F sp h-M5oo for that mass range are 1.84, 1.78 
and 1.73, as compared to 1/qtmYx - 1-78 ignoring the effect. 
The effect is fully negligible for the F sp h(^50o)-A^5oo relation, as 
found above directly from the data; it is at most equal to the sta- 
tistical uncertainty on l/afMY x = 1.78 ± 0.06 (Eq.|2]i and we will 
neglect it in the following. 

In that case, and combining Eq.l2Tl[7land[T6l the Y^-Msoq 
relation for an integration radius of xRsoq can be written as: 



h(zr 2/i Y svh (xR 50 o) = A x 



where 



M 500 



3xl0 14 h^M<» 



a = 1.78; A x = 2.925 x lO' 5 1(x) h^Mpc 2 



I(x) 



= I 3p(u)u z 
Jo 



du 



(22) 

(23) 
(24) 



with p(m) from Eq. [TT| and Q~2] Numerical values for I(x) of 
particular interest are 1(1) = 0.6145 and 1(5) = 1.1037. The 
former gives the normalisation of the F sp h(^50o)-A^500 relation, 
log(A x ) = -4.745. It is in excellent agreement (1% difference) 
with the normalisation derived from a direct fit to the data 
(Eq.[20l). The latter gives the normalisation, log(A x ) = -4.491, 
of the relation for the total Y^ Z D 2 A signal, assuming a cluster ra- 



Similarly, the relation for the SZ signal within an aperture of 
xRsqu is obtained from Eq. [T5l[T3l and[T6l 



h(zr 2 ' 3 Y sz (xR 50 o)D 2 A = B % 
with 



M 500 



3 x lO^h^Mo 



(25) 



a = 1.78; B x = 2.925 x 10~ 5 J(x) Mpc 2 (26) 



J(x) = 7(5) • 



J" 3 p{u) Vm 



2 _ v-2 



x 2 u du 



(27) 



for a cluster extent of 5R^qo. For an aperture of R^q, J (I) - 
0.7398 or log(fi x ) = -4.665. The corresponding Y SZ D A -M 500 
relation is plotted in Fig. [10] We also show measurements for 
A1835 (z = 0.25), A1 914 (z = 0.17), and CL J 1226.9+3332 
(z — 0.89), derived by Mroczko wski et al. (20091) from a joint 
analysis of SZA and X-ray observations using a GNFW pres- 
sure profile model. Although the measurement errors are still 
large, the consistency with the present scaling relation is an en- 
couraging sign of the validity of our determination of the scaling 
relations. Since the clusters cover a wide redshift range, it further 
suggests a standard self-similar evolution, as assumed in Eq. [25] 
Uncertainties on the above relations, that are established 
combining observational and theoretical data, cannot be assessed 
rigourously. Rough estimates of the statistical errors can be de- 
rived by combining the errors on the F sp h(^50o)-Fx and M500- 
Yx relations, with the latter largely dominant. This gives a = 
1.78 ± 0.06 or a - 1.78 + 0.08, further adding quadratically the 
systematic effect of the pressure self-similarity break discussed 
above. The logarithmic error on the normalisation at the pivot is 
+0.024 (±6%). 

6.3.2. Behavior of F sp h(7?) and comparison with the 
isothermal y8-model 

It is instructive to study in more detail the radial dependence of 
Y S ph- Y sp h(R) varies with radius as I(x) (Eq. l24l with p(u) from 
Eq. QT| and [T2l . By construction its normalisation scales with 
mass as Fx- FigurefTTIshows the variation of F sp h(^) with scaled 
integration radius, normalised to Cxsz^x, so that we are effec- 
tively probing F sp h(/?) at fixed mass. 

At large radii, the integrand in I(x) varies as p(u)u 2 du oc uT 2 
for an outer slope /3 ~ 5. As a result, F sp h(^) converges rapidly 
beyond Rsqq and the total SZ signal is not very sensitive to the 
assumption on cluster extent. Assuming a cluster extent of 47?5oo, 
6/?5oo or even lOOT^oo, rather than 57?5oo, changes the total sig- 
nal by only -2%, +1.3% and +4%, respectively. On the other 
hand, the figure shows the dominant contribution of the exter- 
nal regions to Y sp ^: 50% of the contribution to F sp h(/?5oo) comes 
from R > 0.53/?5oo while the region within O.I.R500 and 0.2/?5oo 
contributes by only 2% and 9.5% respectively. This will be even 
more pronounced for the Fsz signal (integration within a cylin- 
drical volume). 

We also plot Ysp h (R) for the GNFW model obtained 
by iNagai et al. (20071) from Chandra data (for the cor- 
rected parameters, [12 .2, 1.3,0.4, 0,9, 5.01, published by 
I Mroczkowski et al. 20091) . It is slightly larger in the centre, as 
expected from the more peaked nature of the scaled Chandra 
profiles (Sec. 13.5b . The agreements very good in the outskirts, 



dial extent of 5R 



500- 



3 Note, however, that (N aeai et al. 20071) assumed a standard self- 
similar mass scaling of the presure profile. The K-M500 relations derived 
from their profiles would differ from ours in terms of slope. 
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Fig. 11. Spherically integrated Compton parameter, y sp h(K), as a func- 
tion of scaled integration radius, normalised to Cxsz^x- Full line: uni- 
versal GNFW scaled pressu re profile (Eq. [TTI and 112b . Dotted line: 
GNFW model obtained by iNagai et al. (20071) for Chandra relaxed 
(T > 5keV) clusters. Dashed-dotted line: isothermal y8-model with 
/3 = 2/3 and a core radius of 0.2Rsoo- 



as it is for the profiles (Fig. [5]), with a slightly higher assymptotic 
value due the slightly smaller value of /3. 

We also compare with the result obtained with an isother- 
mal /3-model, with (5 = 2/3 and a core radius of 0.2/?5 o 
(I Arnaud et al. 20021) . The difference is only 10% at R 500 but 
the model diverges at high radii. This clearly shows that the 
total Ysz signal derived assuming an isothermal /3-model is 
very sensitive to the assumed extent of the cluster. It will also 
be always overestim ated by such a model, as emphasized by 
Hall man et al. (20071) . As an illustration, assuming a cluster ex- 
tent of 2.03 ^5oo, the top-hat virial radius often used in the lit- 
terature, the jS-model gives a total Tsz signal 1 .7 higher than the 
universal pressure profile. 

6.4. The Y-L x relations 

The scaling between the SZ signal and the X-ray luminosity, L x 
is an important relation for comparing X-ray surveys such as the 
ROSAT All Sky Survey and future or on going SZE surveys, such 
as the Planck survey. The luminosity within R^) and in the soft- 
band [0.1-2.4] keV, most relevant for X-ray Surveys , has been 
estimated for REXCESS clusters bv lPratt et al. (2009ah : here we 
used the values both corrected and uncorrected for Malmquist 
bias. 

Figure [12] shows the corresponding F sp h(^5oo)-ix relations. 
We fitted the REXCESS data with a power law: 



fc(zr 2/ %h(*50o) = C 



10 44 h TO 2 ergs s" 1 



h- ( f Mpc 2 



(28) 



The best fitting parameters are given in Table Q] The intrinsic 
scatter around the relation is important, more than 50%, reflect- 
ing the important scatter, at given Y x , of t he soft band lumin osity 
computed without excising the core (see lPratt et al. 20 09a). The 
best fitting relation is consistent with the relation expected from 
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0.1 1.0 10.0 

h(z)- 7/3 L Xi[0J _ 2 . 4] (10 44 h 7 " 2 ergs/s) 

Fig. 12. The F S p h (J?5oo)-£x relations. L x is the [0.1-2.4] keV luminos- 
ity within Rsoo- Full line: power law relation best fitting the REXCESS 
data for Lx corrected for Malmquist bias (black points). Dotted line: 
relation computed by combining the y sp h(-^5oo)-A^5oo relation derived 
from the universal pressure profile (Eq. I22t and the Lx-M^oo relation. 
Dash-dotted Line: best fitting K sp h(^5oo)-^x for uncorrected Lx- 



Table 1. h(zr 2,3 Y sph (R m )-h(z)- 7/3 L x and updated h(z)-~" 3 L x - 
M500 relations (see text). L x is the [0.1-2.4] keV luminosity 
within Rsoo- MB: relations corrected for Malmquist bias. For 
each observable set, (B,A), we fitted a power law relation of 
the form B = C(A/A ) a , with A = 10 44 h TO 2 ergs/s and 3 x 
10 14 h 7( J M for L x and M500, respectively. cn og y. intrinsic scatter 
about the best fitting relation in the log-log plane. 



Relation 


log 


c 




a 


0"lo 




Y sph (R 50O hL x -MB 


-4.940 


± 0.036 


1.07 


±0.08 


0.190 ± 


0.025 


L X -M 500 -MB 


0.193: 


t 0.034 


1.76 


±0.13 


0.199 ± 


0.035 


^sph(^50o)-ix 


-5.047 


± 0.037 


1.14 


±0.08 


0.184 ± 


0.024 


Lx-M 5 oq 


0.274 = 


t 0.032 


1.64 


±0.12 


0.183 ± 


0.032 



combining the 5 / sp h(/?5oo)-^500 relation derived from the univer- 
sal pressure profile (Eq. l22l and the Lx-A^soo relation. For con- 
sistency, the latter was updated (parameters given in Table [TJ, 
using present M500 values derived from the updated M500-TX re- 
lation (Eq. |2j. The slope and normalisation (taking into account 
the different pivot used) are consistent with those published in 
iPratt et al (2009al) . 

For practical purposes, the scaling of Ysz(xRsoo) D\ or that 
of the total SZ signal with L x is of more direct interest than 
the T sp h(^5oo)-^x relation. In view of the good agreement of the 
latter with the universal profile model, the Y-L x relation, for 
any integration region of interest, can be safely derived by cor- 
recting the normalisation in Eq. [28] by the model ratio of Y to 
T sp h(^50o)- This ratio is simply I(x)/I(l) for the spherically inte- 
grated Compton parameter, e.g., 7(5)//(l) = 1.796 for the total 
SZ signal, and J{x)/I{l) for the Tsz(*^50o) D 2 A signal. 
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6.5. Comparison with standard self-similar relations 

The F-M500 relations derived above do not seem to deviate much 
from standard self-similarity (Fig. ITOb . A fully consistent stan- 
dard (ST) model, with standard slope T-M500 relations, is ob- 
tained when using the standard slope A/500-Tx relation (Eq. [3]), 
as shown in Appendix IE1 The universal profile and scaling rela- 
tions obtained in that case are given in the Appendix, together 
with a detailed comparison of the presently derived scaling re- 
lations with the ST relations. In summary, the difference for the 
F-M500 relations mirrors that for the M^qq-Yx relation. As com- 
pared to values derived from the ST relation, Y is lower at low 
mass and higher at high mass. Typically, the difference for the to- 
tal Y sz signal ranges from -19% to +6% in the [10 14 -10 15 ] M 
mass range. On the other hand, the Y-Lx relations, which only 
depend on cluster internal structure, are essentially the same 
in the two models : the difference is less than 5% in the [0.1- 
10] 10 44 ergs/s luminosity range. 



7. Discussion and conclusions 

The present work is the first examination of the properties of 
the ICM pressure for a representative sample of nearby clusters. 
The sample, REXCESS, was chosen by X-ray luminosity alone, 
without regard to morphology or dynamical state. It covers the 
mass range 10 14 < M500 < 10 15 M , with mass iteratively esti- 
mated from the M500-TX relation, calibrated from a sample of 
relaxed clusters in cluding REXCESS objects. As for the entropy 
dPratt et al. 2009bl) . the depth of the observations allowed us to 
probe the scaling behavior of the pressure profiles out to R^qq. 
This is essential for a complete picture of the modification of 
the standard self-similarity due to non-gravitational processes, 
including its radial behavior. 

Scaling the individual pressure profiles by mass and redshift 
according to the standard self-similar model, we derived an av- 
erage scaled pressure profile for the cluster population and relate 
the deviations about the mean to both the mass and the thermo- 
dynamical state of the cluster: 

- Cool core systems exhibit more peaked profiles, while mor- 
phologically disturbed systems have shallower profiles. 

- As a result, the dispersion is large in the core region, reach- 
ing approximately 80 per cent at 0.03 Rsqq. However, as com- 
pared to the density, the pressure exhibits less scatter, a result 
of the anticorrelation of the density and temperature profiles 
interior to 0.2/?soo. Outside the core regions, the dispersion 
about the average profile is remarkably low, at less than 30 
per cent beyond 0.2 R500. 

- We find a residual mass dependence of the scaled profiles, 
with a slope of ~ 0.12, consistent with that expected from 
the empirical non-standard slope of the M^-Yx relation. 
However, there is some evidence that the departure from 
standard scaling decreases with radius and is consistent with 
zero at R50Q. We provide an analytical correction to the mean 
slope that accounts for this second order effect. 

The behaviour of the pressure profiles, with respect to stan- 
dard self-similarity with zero dispersion, resembles that gener- 
ally found for other quantities such as the entropy or density: 1) 
regularity in shape outside the core 2) increased dispersion in- 
side the core linked to cooling effects and dynamical state and 
3) departure from standard mass scaling that becomes less pro- 
nounced towards the cluster outskirts. However, the latter two 



deviations are less pronounced than for the entropy and/or den- 
sity, showing that the pressure is the quantity least affected by 
dynamical history and non-gravitational physics. This further 
supports the view that isz is indeed a good mass-proxy. 

Furthermore, our direct measure of the T sp h(^5oo)-Tx 
relation, where T sp h(^50o) is the spherically integrated pressure 
profile, exhibits dispersion consistent with the < 5% statistical 
scatter. This shows that variations in pressure profile shape 
do not introduce significant extra intrinsic scatter into the 
T S ph(^50o)-A^500 relation as compared to that from the Tx-^500 
relation. 

The observational data are compared to and combined with 
simulated data to derive the universal ICM pressure profile. 
This profile is then used to predict the scaling relations in- 
volving the integrated Compton parameter Y. We consider both 
the spherically integrated quantity, F sp h(/?), which is related to 
the gas thermal energy, and also the cylindrically integrated 
quantity Y cy i(R) = Ysz(R)D 2 A , which is directly related to the 
observationally-derived SZ signal within 6 = R/Da '. 

- Simulated scaled profiles from three independent sets of state 
of the art numerical simulations show excellent agreement, 
within 20%, between 0.1 and 3/?5oo, for pressures varying by 
4 orders of magnitude in that radial range. 

- Comparison with observed scaled data shows good agree- 
ment outside the core regions, which is the most relevant as- 
pect for the Fsz estimate. The average simulation profile lies 
parallel to the observed data, with only a slight offset (~ 10 
per cent) when the simulated profiles are scaled using the 
hydrostatic mass. 

- This motivates us to combine the average observed scaled 
profile in the [0.03 - l]/?5oo radial range with the average 
simulated profile in the [1-4] R500 range. This hybrid pro- 
file is fitted by a generalised NFW model, which allows us 
to define a dimensionless universal ICM pressure profile. 
Combined with the empirical mass scaling of the profiles, 
this universal profile defines the physical pressure profile of 
clusters, up to the cluster boundary, as a function of mass and 
redshift, assuming self-similar evolution. 

- This universal profile allows us to derive the expected 
T sp h(x#5oo)-M 5 oo or F sz (x# 5 oo)-^500 relations for any aper- 
ture. The slope is the inverse of the empirical slope of the 
Msoo—Yx relation. The normalisation is given by the dimen- 
sionless integral of the universal profile within the region of 
interest expressed in scaled radius. The corresponding Fsz- 
Lx relations can be derived by combining the relevant Y— 
M500 relation with the empirical Lx- M500 relation. 

- The Y sph (R 500 )-M 500 and 5 / sp h(^ 5 oo)- L x relations derived di- 
rectly from the individual profiles are in excellent agreement 
with those expected from the universal profile. 

- We confirm that the isothermal yS-model over-estimates the 
Y signal at given mass. This overestimate depends strongly 
on the assumption on cluster extent and reaches a factor of 
nearly two at 2/?5oo- 

The convergence of various approaches to determine scaled clus- 
ter profiles supports the robustness of our determination of the 
universal pressure profile, particularly its shape. This includes 
the agreement between independent simulations, between these 
simulations and the present observed data based on a representa- 
tive cluster sample, and also the agreement between the present 
XMM-Newton data and published Chandra data for clusters of 
similar thermo-dynamical state. As a result, we believe that 
quantities which purely depend on the universal profile shape are 
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particularly robust and well converged. This includes the typical 
SZ decrement profile or relations between the Compton param- 
eter estimated in various apertures. 

However, the pressure profile interior to Rsoq is derived 
from temperatures estimated using azimuthally averaged spec- 
tra. These have been corrected for the spectroscopic bias due to 
projection but not for azimuthal variations. In the cluster out- 
skirts, the el ectron-proton eq uilibration time is larger than the 
Hubble time dFox & Loeb 1997b and if the electron temperature 
is indeed smaller than the ion temperature, this will affect the 
pressure profile and l ead to a decrease in the total Y sz signal 
dRudd & Nagai 20091) . High resolution SZ data with improved 
sensitivity are needed to probe any remaining systematic effects 
due to the spectroscopic bias, and to directly observe the shape 
of the pressure profile beyond Rsm, which is out of reach of cur- 
rent X-ray observatories. 

Using the universal profile, the absolute normalisation and 
slope of the Y-M500 relations rely on the underlying obser- 
vationally defined M^oq-Yx relation. Initial comparison with 
FszCffsoo) data for 3 high m ass systems, measured with SZA by 
Mro czkowski et al. (20091) and analysed with a realistic analytic 
pressure profile, indicates good agreement. A key point is to ex- 
tend this type of analysis to larger samples and to include lower 
mass systems. We further emphasize that the Msorj-Ix relation 
was calibrated from hydrostatic mass estimates using relaxed ob- 
jects. The Y-M500 relation we derive is technically a Y-X-ray 
mass relation and is expected to differ from the 'true' F-M500 by 
the offset between the 'true' mass and the hydrostatic mass for 
relaxed objects. 

A major open issue is the pressure evolution. With the 
present study based on a local cluster sample, we could only 
assume standard self-similar evolution. Because the SZ signal 
is not subject to redshift dimming, on going SZ surveys are 
expected to detect many new clusters at high z. Of particular 
interest is the Planck survey, which, thanks to its All-Sky cov- 
erage, will detect massive, thus rare, clusters, the best objects 
for precise cosmology with clusters. SZ follow-up, at the best 
possible resolution, and sensitive X-ray follow-up (particularly 
with XMM-Newtori) will be crucial to assess possible evolution 
of pressure profile shape and measure the evolution of the 
M^-Yx and Fsz-^500 relations. 

As a matter of practical application, the universal pressure 
profile is given in Eq.QT|with parameters in Eq.[T2] For clusters 
of given mass M500 and z, the physical pressure profile can then 
be derived from Eq. [13] and the spherical Y sp h{R) or cylindrical 
Y cy \(R) quantities can be estimated for any radius of interest us- 
ing Eq. l22H24l and Eq. [2514271 respectively. 

These equations can be used as is when M500 is estimated 
for relaxed systems using the HSE equation, and for all clusters 
using M500 derived from mass-proxy relations. The preferential 
relations would be the Msm-Yx and the M sqq-Lx, where L x is 
the core-excised bolometric luminosity dPratt et al. 2009al) . as 
both these relations display low scatter, compared to the rela- 
tion between M500 and the full aperture soft band Lx- A typical 
application would be to predict the SZ signal of a known X-ray 
cluster with measured Lx or M500, or to estimate the mass and 
thus X-ray properties of newly discovered SZ clusters. Other ap- 
plications include the analysis of low S/N and/or poor resolution 
SZ observation of X-ray clusters, e.g., allowing to optimise the 
integration aperture and use a realistic decrement shape. 

On the other hand, care is needed when knowledge of the 
'true' mass is important, e.g., in predicting cluster number 
counts for future SZ surveys or in SZ selection function mod- 



elling. The above total Fsz-A^oo relation should be corrected by 
the bias between the true mass and the HSE mass at R500, which 
is typically ~ 13% as determined from comparison with cur- 
rent numerical simulations. Further progress on this fundamental 
question, as well as on the intrinsic scatter of the Y-M relation, is 
expected from the wealth of high quality multi-wavelength data 
that will be available in the coming years. 
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Appendix A: Characteristic self-similar quantities 

Following IKfagai et al. (20071) and I Voit et al. (2005b the charac- 
teristic quantities, P500 and Y500, used in the present work, are 
defined from a simple self-similar model. The characteristic 
temperature is k^oo = nm v G M^/lRsca, the temperature of 
a singular isothermal sphere with mass M500. Here, ft is the 
mean molecular weight and m p , the proton mass. We recall that 
M500 is defined as the mass within the radius R^) at which the 
mean mass density is 500 times the critical density, p c (z), of the 
universe at the cluster redshift: M500 = (An / '3) R^ 00 500 p c (z) 
with p c (z) = 3//(z) 2 /(8ttG). H(z) is the Hubble constant, 
H(z) = H(0) VQm(1 +z) 3 + Q a and G is the Newtonian con- 
stant of gravitation. The characteristic gas density is p g ,5oo = 
500/b/5 c (z), i.e., the ratio of the gas density to the dark matter 
density is that of the Universe baryon fraction /b . The electron 
density is ra e ,50o = Pg,soo/(He m?) where ju e is the mean molecular 
weight per free electron. 

The characteristic pressure, P500, is then defined as: 



^500 - «e,500k7'500 

" 500 G~ 1/4 H(z) 2 14/3 



3_ 

8^ 



^ f M 2/3 



(A.l) 
(A.2) 



and the corresponding characteristic integrated Compton param- 
eter is: 



T500 



(Tj 

m e c 2 



— R 5m P500 



cr T f B M500 kr 5 oo 



/500GH(z) 



m e cr 
2/3 



He 



(A3) 
(A.4) 



Numerical coefficients given in the corresponding Eq.[5]and 
[16] are obtaine d for f H = 0.175, h = 0.59 and fi e = 1.14, the val- 
ues adopted by Nag ai et al. (20 07). allowing a direct comparison 
with their best fitting GNFW model. Note that the exact choice 
for these parameters does not matter, and does not need to reflect 
'true' values, as long as the same convention is used throughout 
the study (e.g., when comparing observed and theoretical scaled 
profiles or observed scaled profiles from different samples or in- 
struments). 
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Fig. B.l. Ratio of the scaling relations derived using the empirical 
M 500 -Y x relation (Eq.[2} to those derived using the standard slope rela- 
tion (Eq.O. From top to bottom: M500 as a function of Yx', Y\, 1^11(^500) 
and total Y sz as a function of M500; L x as a function of M 500 ; KsphC^soo) 
and total Ysz as a function of Lx. 



words, the pressure profiles follow a standard self-similar mass 
scaling: 



P(r) = P 5 oop(r/R 5m ) 



(B.l) 



Appendix B: The standard self-similar case 

In this Appendix, we summarise results (hereafter ST results) 
obtained when M500 is estimated for each REXCESS clusters 
using the Msoo-ix relation with a standard slope (Eq. [3}. The 
other physical parameters are consistently estimated, R500, Yx 
and Tx simultaneously in the iteration process used to derive 
M500 (Sec. and T S ph(^50o) from integration of the pressure 
profiles up to ^500- For practical purposes, the baseline param- 
eters obtained using the best fitting empirical Msoo-^x relation 
(Eq.[2ji can be converted to the ST values using the power law re- 
lations given in Table IB. ll The luminosity Lx is kept unchanged, 
the difference in R500 values (at most 4.5%), having a negligible 
impact due to the steep drop of emission with radius. 

In the ST case, the scaled pressure profiles do not show any 
significant dependence on mass, as shown in Sec. 13.41 In other 



Table B.l. Power law relations to convert physical parame- 
ters of REXCESS clusters from those derived using the empir- 
ical Mgoo-Tx relation (Eq. |2]i to those derived using the stan- 
dard slope relation (Eq.[3]l. For each observable, Q, the conver- 
sion follows the form Q ST = C(Q/Qo) a where the pivot, Qo 
is 3 x 1O 14 M , 5keV, 2 x lO 14 M keV and 2 x lO^Mpc 2 for 
M 5m ,T x , Y x and T sph (^50o), respectively. 



Relation 


C 


a 


Y sph (Rsoo)~Y S p h (R 5QO ) 


0.968 
1.002 
0.995 
0.991 


1.089 
0.992 
1.017 
1.031 
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Table B.2. h(z)- 2/3 Y s?h (R 500 )-h(zT in L x and h(z)- 7/3 L x -M 50 o 
relations for M500 estimated using the standard slope M500-FX 
relation (Eq.0). Same notations as in TableQ] 



Relation 


log 


10 c 


a 




trio 




F sph (i?5oo)-ix-MB 


-4.947 


± 0.037 


1.08 ± 


0.08 


0.192 ± 


0.025 


L X -M 5Q0 -MB 


0.215 = 


t 0.035 


1.61 ± 


0.12 


0.199 ± 


0.035 


?sph(^50o)-ix 


-5.056 


± 0.038 


1.16 ± 


0.08 


0.184 ± 


0.024 


Lx-M 5Q o 


0.295 - 


t 0.032 


1.50 ± 


0.11 


0.183 ± 


0.032 



Appendix C: Pressure profiles and best fitting 
model 

Here we list the physical cluster properties and the parameters of 
the GNFW model best fitting each profile (Table [Ct. Individual 
profiles and their best fitting model are plotted in Fig.Q~|-[3] 



with P500 defined by Eq. [5] The GNFW parameters of the uni- 
versal profile p(x), derived as described in Sec. [5] are: 



[P ,c 500 ,y,a,p] = [8.130h TO 3/2 , 1.156,0.3292, 1.0620,5.4807] (B.2) 



As a result, the integrated Compton parameters also follow 
standard self-similarity, F oc Ml£. The F-M500 relations derived 
from the universal pressure profile can be written as: 

Y svb (xR 50 o) = Y 500 I(x) (B.3) 
Y sz (xR 500 )D 2 A = Y 50Q J(x) 

with F500 given by Eq. [16] and I(x) or J(x) defined by Eq. [24] 
and Eq. [27] respectively. For the GNFW parameters given by 
Eq. IB. 21 the numerical values of 1(1), 1(5) and 7(1) are 0.6552, 
1.1885 and 0.7913, respectively. The F sp h(x/?5oo)-Mg o relation 
derived from a direct fit to the data has a slope of 1.663 ± 0.013, 
fully consistent with 5/3. Over the [10 14 — 10 15 ] M Q mass range, 
it differs by less than 0.8% from that derived from the universal 
profile (Eq.|B3T>. 

We also derived the observed T S ph(^50o)-Lx relation, as well 
as the Lx-M soq corresponding to the modified M500 values. The 
best fitting power law parameters are given in Table IB. 21 The 
former is consistent with the relation expected from combining 
the Lx-^500 relation with the T sp h(^5oo)-M5oo relation derived 
from the universal pressure profile (Eq. IB.3b . The Y-Lx relation, 
for any integration region of interest, can be derived by correct- 
ing the normalisation of the F sp h(x/?5oo)-£x given in Table IB~2l 
by the model ratio of Y to F sp h(/?5oo), as described in Sec. 16.41 

Figure iBTI compares the scaling relations derived in the pa- 
per with the ST relations derived in this section. The empirical 
slope of the M^qq-Yx relation being smaller than the standard 
value, M500 at a given Fx is higher at low Yx and smaller at high 
Fx (top panel). Equivalently, Fx at given mass is smaller at low 

mass, by 16% at M500 = 10 I4 M o , and higher at high mass, 

by ~ +10% at M500 = 1O 15 M (second panel). The behavior of 
Fsz closely follows that of Fx (same panel) simply because the 
ratio of the two purely depends on the shape of the universal pro- 
file. This shape is barely affected by the small difference in Z?soo 
values used to scale the physical pressure profiles. Similarly, the 
Fsz-^x relation only depends on cluster internal structure and is 
essentially the same in the two models (bottom panel). Fsz(^50o) 
is slightly higher/lower at low/high Lx following the change of 
/?5oo at given Lx- As the Mgoo-Fx is shallower than the ST rela- 
tion, the A/500-ix is also shallower (thus higher R500 at low Lx) 
or equivalently the LX-M500 is steeper (third panel). 
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Table C.l. Cluster physical parameters. Column (2)-(3): R^m is the radius corresponding to a density contrast of 500, estimated 
iteratively from the Msqo-Ix relation (Eq.|2]i, where Yx = M g; 5oo7x is the product of the gas mass within R500 and the spectroscopic 
temperature Tx- Column (3) and (4): spherically integrated Compton parameter within R2500 and R500, respectively. Column (5): 
P500 as defined by Eq. [5] Column (6) to (9) give the best fitting GNFW parameters for the pressure profi les (Eq.[TT[). The e xternal 
slope parameter/? has been fixed to 5.49 (see text). Redshift z and M500 values can be found in Table 1 of lPratt et al. ("2009bh . 



Cluster 




Y x 


^sph(^250o) 


^sph 




P 5QQ 


^0 
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3 keV) 


(10 


MpC ) 


O0~ 5 
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Fig. 2. continued 
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Fig. 3. continued 



